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We extend a theorem of Bloch, which concerns the net orbital current carried by an interacting 
electron system in equilibrium, to include mesoscopic effects. We obtain a rigorous upper bound to 
the allowed ground-state current in a ring or disc, for an interacting electron system in the presence 
of static but otherwise arbitrary electric and magnetic fields. We also investigate the effects of spin- 
orbit and current-current interactions on the upper bound. Current-current interactions, caused by 
the magnetic field produced at a point r by a moving electron at r', are found to reduce the upper 
bound by an amount that is determined by the self-inductance of the system. A solvable model of 
■ an electron system that includes current-current interactions is shown to realize our upper bound, 

0^ ' and the upper bound is compared with measurements of the persistent current in a single ring. 
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I. INTRODUCTION 

A theorem due to Bloch holds that an interacting electron system in equilibrium carries no net orbital current 
t^J- ■ [0. This question originally had been motivated by early attempts, before the Bardeen-Cooper-Schrieffer theory, to 
explain superconductivity by proposing that electron-electron interactions lead to special current-carrying states of 
lower energy than the current-free states. However, it is now understood that supercurrent-carrying states are in fact 
metastable nonequilibrium states, which, because of their off-diagonal long-range order or wave function rigidity, have 
an extremely long lifetime. 

The great interest over the past several years in the physics of mesoscopic systems again has made the question of 
allowed equilibrium currents an important one. More than ten years ago, Buttiker et al. Q predicted the existence 
i of equilibrium currents in mesoscopic normal-metal rings threaded by a magnetic flux. Recent experimental evidence 
HQ in support of this conjecture has stimulated considerable interest in these so-called persistent currents. Although 
a satisfactory explaination of the experiments is still lacking, the present consensus is that both electron-electron inter- 
action and disorder effects are important g . A related phenomena is that of spontaneous orbital currents occurring in 
the absence of any applied magnetic field or twisted boundary conditions. Although there is no experimental evidence 
for this symmetry-breaking state, spontaneous orbital currents have been predicted to occur by several authors 

Given the diverse situations in which equilibrium current-carrying states may occur, it is worthwhile to revise 
Bloch's theorem to incorporate these mesoscopic effects. To this end, Vignale || has recently derived a rigorous 
upper bound to the persistent current in a single ring, the results being valid for both noninteracting and interacting 
electrons in the presence of arbitrary magnetic fields and impurity potentials. One surprising result of Vignale's 
analysis is that although the upper bound to the persistent current in a thin ring of uniform density and radius R 
vanishes as 1/R for large R, it does not vanish for a thick ring or punctured disc with the ratio Ri n /Rout of the inner 
radius and outer radius fixed as these radii become infinite. Although Vignale's result does not preclude the existence 
of a more stringent upper bound that always vanishes in the macroscopic limit, the upper bound is actually realized in 
calculations of the persistent current (the integrated azimuthal current density) in a two-dimensional noninteracting 
electron gas in a large quantum dot pfj| , pT| . 

One motivation for this work is to extend the analysis of Ref. || to include the effects of spin-orbit interaction, 
which has received considerable attention in connection with persistent currents and spontaneous currents. Spin-orbit 
interaction is known to lead to a topological interference effect, called the Aharanov-Casher effect [j^], which is an 
electromagnetic dual of the Aharanov-Bohm effect. Meir et al. Jll| have shown that spin-orbit scattering in one- 
dimensional disordered rings induces an effective magnetic flux which reduces the persistent current in a universal 
manner. The effect of spin-orbit interaction on mesoscopic persistent currents has been studied by several other 
authors, who also find reduced currents This has led us to question whether the upper bound on the allowed 

persistent current is itself reduced by spin-orbit interactions. We shall show here that this is not the case. 

A second motivation for this work is to examine the infleuence of current-current interactions, an order v 2 /c 2 
relativistic effect caused by the magnetic field produced at a point r by a moving electron at r'. The possibility 
of these magnetic interactions leading to a spontaneous current-carrying state in mesoscopic metal rings has been 
discussed in a remarkable paper by Wohlleben et al. p], where it is shown that, in zero field, a small ring exhibits 
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a transition to a state with persistent current. The combined effects of spin-orbit coupling and current-current 
interactions has been studied recently by Choi We shall show below that current-current interactions reduce the 
upper bound on the ground-state current by an amount which is determined by the self-inductance of the system. 

In this paper, we derive an upper bound to the ground-state current in an arbitrary many-electron system, including 
spin-orbit coupling and current-current interaction effects. To best demonstrate the modifications to Bloch's theorem 
from the effects of finite sample size, we restrict our analysis to zero temperature. However, our final results are also 
valid at finite temperature, as may be shown by following the method of Ref. 0. 



II. RIGOROUS UPPER BOUND 



We begin by obtaining a many-electron Hamiltonian that includes current-current interactions to order v 2 /c 2 . In 
the transverse gauge, the vector potential seen by an electron at r n in the presence of the other moving electrons (of 
charge — e) is 

where v n is the velocity of the nth electron, and where 

T^{r) = S ij + (2) 
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This vector potential leads to a current-current interaction term in the Lagrangian of the form 

2c 2 



e 2 ^ vj r' 3 (r n — iy) iP n , 

L ^ -^2^ lr_ -r_,l ' (6) 



n<n' 



which in turn leads to a current-current interaction term in the Hamiltonian of the form 



n<n' 

to leading order in v 2 /c 2 . The complete Hamiltonian, including spin-orbit coupling, Coulomb and current-current 
interactions, and coupling to additional electric and magnetic fields, may be written as 

" = E (S - + "<'•> + ' l w(r »» ' n »0 + E k^I - is? E n »n^^r n - 

n x ' n<n' 1 1 n<n' ' ' 

(5) 

where n„ = p„ + -A(r n ), the S* 1 are spin operators, and where the potentials A and V are time-independent but 
otherwise arbitrary. Spin-spin interactions and the coupling of the spin degrees of freedom to the external magnetic 
field are not important here and shall be ignored. The velocity operator v„ = [r„, H]/ih is given by 

t W n ( _Hj_\ , 1 ro w t e 2 v - / r^(r»-iy) T' J " (r» - r„- ) ■ \ 

We shall consider a system of iV electrons confined to a ring or disc, oriented with its axis along the z direction, and 
we write (0) in cylindrical coordinates r = (r, 8, z). The thickness and cross-sectional shape of the system is arbitrary. 
The many-body ground state ip(r 1 s 1 , ■ ■ ■ ,r N s N ) satisfies 

Hip = Eil), (?) 

where E is the ground-state energy 

Suppose that the ground state ip carries an orbital persistent current 



E 



ee(r„) e g (r n ) 
v„ + v„ 



1>), (8) 
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where eg(r) is an azimuthal unit vector at r. We may construct a rotating state ip' = Utp, where 

U = Y[e lSLe " /h , (9) 



which is not necessarily an eigenstate of H, and which has a mean energy E' = (ip'\H\ip') given by 

->2 



1 n} 1 



4m 2 c 2 \ ' r n |r„ — r„/ |r„/ / 2m 3 c 2 \ ^ r 3 / 8m 3 c 2 

Here "( )" denotes an expectation value in the original ground state ip, H 9 n = II„ • eg(r n ), and we have used 

£/ t n„£/ = n„ + — e e (r n )6L. (11) 

The energy difference, SE = E' — E, plotted as a function of the parameter SL, is shown in Fig. 1. For values of 5L 
such that 

< SL < SL*, (12) 

where SL* is the zero of 8E defined in Fig. 1, the rotating state ip' has a lower mean energy than ip, so ip cannot be 
the ground state. 

This is the essential content of Bloch's theorem. It applies whenever there is a nonzero SL satisfying ([l2]). However, 
the smallest nonzero SL permitted by the condition that the wave function ip' be single-valued is SL = h. When the 
relativistic corrections in (|E|) are neglected, SL* = SLq, where 

Airml I 1 \ 1 

dL = 

Here 



i\ s !/v^\=l/>r=« (14) 
ri / rl N J r"i ' y ' 

where n(r) is the ground-state number density. Therefore, Bloch's theorem applies only when SL* > h, or whenever 
1^1 > J£ax, wh ere 



This is the upper bound derived in Ref. ||. When the relativistic corrections are included to leading order, SL* is 
given by 

SL = SL* + emc 2 N 2Qi r 2}2 ( ( n « 71 + 72" n » + 4 ~^T~ 



n n 



2lTeI I 4( r „) r^(r»-ry) e'(iy) \ 16tt 2 I 2 / v Dg \ 167rW(l/r 4 ) 

c 2 iV 2 (l/r 2 ) 2 \^ ( r„|r n -r n ,|r n , / e 2 c 2 iV 3 (l/V 2 ) 3 \ ^ r 3 / e 3 c 2 iV 3 (l/r 2 ) 4 ' 1 ' 

Bloch's theorem therefore applies whenever |/| > / max , where 

/max = (1 - A)/° ax , (17) 

and 
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1 /yiV 1 | 1tt 2 i j^HN i p2 / V 4(* n )TV(r n -T n ,)e? e (r nl y 

4m 2 c 2 N(l/r 2 ) \^\ n r% r% n r\ ) / 2mc 2 N(l/r 2 \ \ ^ r_lr_ - r_.lr_. 



m 2 c 2 ^(l/r 2 ) \ ^ r,3 / 4m 2 c 2 (l/r 2 ) V ' 

is a dimensionless reduction factor. States carrying orbital currents larger than / max cannot be ground states of (^). 

The upper bound (|l7| ) applies to interacting electron systems in the presence of static but otherwise arbitrary 
electric and magnetic fields, and includes the effects of spin-orbit coupling and current-current interaction. The upper 
bound (|l5|) applies to noninteracting systems and also to electrons with Coulomb interaction. In particular, ( |15| ) 
applies to noninteracting electrons in a periodic potential, and this fact leads to a general constraint on the band 
structure of any one-dimensional crystal |lq] . 



III. UPPER BOUND FOR A THIN RING 



Now consider the case of a thin ring with cross-sectional dimensions much less than the radius R of the ring. In 
this case, 

Nek 2e^ 

raax AirmR 2 L ' [ ' 

where v F is the Fermi velocity, L = 2nR is the circumference of the ring, and 

A* 3 /V(tff\ + A/ y e|(r n )rO-(r n -r n 4(r B Q \ h /y n e\ + g (2Q) 

2m 2 c 2 JV\^ 1 n > I 2mc 2 N\^ |r n -r„/| / m 2 c 2 NR\^ " / 4mW V ; 

The first term in (|2^) is approximately equal to E F /mc 2 , where E F is the Fermi energy, and hence this term is 
entirely negligible here. The magnitude of the third term may be estimated by using the approximation ^ ^2 n H^) ~ 
47rmi?/^ lax /e, which shows that the third term and fourth term in ( |20| ) are both of order X 2 /R 2 , where A c = h/mc is 
the Compton wavelength of the electron. These terms are therefore negligible here as well. 
The operator in the second term of ( po| ) may be written in second-quantized form as 

^ e» 9 (r n ) T*(r n -r „,) eP e {v n ,) = f^^p^^^^^ (21) 
where tp(r) and ip ( r ) are electron field operators, and where 

FMs mi^iAn. ( 22) 

This term is a consequence of the current-current interactions. In a mesoscopic or macroscopic system, the largest 
contribution to the expectation value of (|l]) comes from the direct term 

d 3 rd 3 r' F{r,r') n(r)n(r'), (23) 

which is normally absent in the case of Coulomb interactions in a uniform system. For a thin wire of approximately 
uniform density and current density, we may write the latter as I divided by the cross-sectional area N/nL, 

j(r) » ^ e e (r). (24) 



Then we have 



\- ^(r»)T^(r„-r»Qe J 9 (r„0 \ _ iV 2 f 3 3 , /(r)T^(r - r')f(r') 

2^ lr_ -r_,l / - 4^2 R 2r2 / d rCt r 



, 4tt 2 R 2 I 2 J |r-r' 
N2 /,3„^jW-jt(r') 



tfV dV (25) 



2vr 2 i? 2 / 2 / |r-r 
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where jt is the transverse current density, defined as 



j t (r) = x V x J d 3 



|r — r' 

i „ {*,v>m 



:j(r) + — V / rfV J \ n < . (26) 



Because the equilibrium current density satisfies V • j = 0, the second term in (]26J) vanishes and j is purely transverse. 
The reduction factor (|2(]) for a thin ring may therefore be written approximately as A = 2/" lax £/c$0i where 

c ^ i [ d *r d \> £tm (27) 

I 1 J |r-r'| 

is the classical self-inductance of the ring, and where $o = hc/e is the quantum of flux. It is also useful to rewrite the 
reduction factor as 

21° 

A = -iss, (28) 

where I c = c$o / C \s the magnitude of the current needed to produce one quantum of flux. This latter form makes 
explicit the relative importance of the inductive effects. Therefore, the upper bound in a thin ring may be written as 

27"° \ 

(29) 

We see that the current-current interactions always reduce the allowed ground-state current by an amount that depends 
on the self-inductance of the ring. This reduction occurs because the energy required to sustain a persistent current 
now includes the magnetic field energy. As is clear from our derivation, which treated the current-current interaction 
as a small perturbation, ( |29| ) is valid only when A << 1 , or when I^ ax « I c . We shall evaluate (p9| ) for realistic 
thin-ring geometries in the final section of this paper. 

IV. PERSISTENT CURRENT IN A SOLVABLE MODEL WITH CURRENT-CURRENT INTERACTIONS 

In this section we shall calculate the ground-state current in an electron gas with current-current interactions, which 
is confined to a thin wire loop of circumference L = 2itR. We shall assume that only a single transverse mode in 
the wire is occupied, so that the system is effectively one-dimensional, with a width a « L. For modest applied 
perpendicular magnetic fields, such that the magnetic length I = (hc/eB)? satisfies £ » a, the primary effect of 
the magnetic field on the electron gas is to induce an Aharanov-Bohm phase shift around the ring. The effect of the 
magnetic field in this regime may therefore be accounted for by imposing the twisted boundary conditions 

Mx + L) = e^Un(x), (30) 

on the single-particle states. Here x = R9 is an arclength coordinate going around the circumference of the ring, and 
4> = $/$o is the total enclosed flux in units of the flux quantum <I>o = hc/e. 

We shall calculate the ground-state current in a self-consistent mean field approximation in which the effect of 
the current-current interaction is to increase the flux enclosed by the ring by an amount that is determined by the 
persistent current. The persistent current is, in turn, determined by the total flux. Coulomb interactions, which are 
not expected to affect the persistent current in a disorder-free system, are ignored altogether. The single-particle 
states are therefore of the form 

^n = ^= e lk - x , (31) 

where n is an integer and, according to (|30|), the allowed wave numbers are k n = 2ix(n + <fi)/L. The energies of the 
states ( |3l| ) are e„ = h 2 k^/2m, and 

/ n = -^(n + 0) (32) 
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are their contributions to the total persistent current. 

Let N be the total number of electrons, which we shall take to be twice an odd integer. (A similar analysis applies 
to the case of N/2 even.) The total persistent current at zero temperature is given by 

I = 2^(E F - e n ) I n , (33) 

n 

where E F is the Fermi energy, 9{x) is the unit step function, and where the factor of 2 accounts for the spin degeneracy. 
For a small flux, the set of occupied states in ([33|) is not changed by the flux. Then we find 

/ = -2/° ax 0, (34) 

which is the well-known result for the persistent current in a one-dimensional ring in the presence of an external 
dimensionless flux 0, when disorder and electron-electron interactions are ignored. The expression (|34|) is valid until 
the set of occupied states is changed by the flux. In the usual case where cf> is independent of I, the range of allowed 
flux is — | < <j> < i , and the upper bound ([n]) is realized when — > ± i . 
In our case of interest, however, there are two contributions to (f>, 

= 01 + <t>2- (35) 

The first, <pi, is the flux from the external magnetic field, and the second, 4>2, is the flux that originates from the 
current-current interactions. The latter is the flux through the ring induced classically by the current /, 

$ 2 = - / rfa-Vx / d 3 r'j^-. (36) 
cJring J |r-r'| 

In a thin ring, this leads to 02 = IC/c&o, where £ is the self-inductance. Therefore, we obtain from (|34|), 

7 = -2/° ax 0i(l-2%^). (-i<0 1 <|) (37) 

In this simple model, the current is zero when the external flux is zero. The maximum current occurs when </>i — > 
The magnitude of the persistent current in this optimal case is then found to be 

\I\=^Jl-^^), (38) 



which realizes our thin-ring upper bound (E9J) for electrons with current-current interactions 



V. DISCUSSION 

We now evaluate the upper bound (|2^) for realistic thin-ring geometries. The upper bound ( |l9| ) for a metal with a 
Fermi velocity of 2x 10 8 cm/s, a typical value, may be written as 

7° ~, °- 64/iA (39) 

-'max t I \ ' V°^> 

L(/jm) 

where L(/Ltm) is the circumference of the ring in microns. As discussed above, the relevance of inductive effects are 
characterized by the ratio of I c = c$o/ which is the current needed to produce one quantum of flux, to 7 max . If we 
measure the self-inductance in microns, then 

41.2 mA . , 

The reduction factor ( |28| ) for a metal ring may then be written as 

A«3.0xI(T 5 y. (41) 

1j 
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The self-inductance of a thin toroidal ring with major radius R and minor radius a (wire radius) is C — 
4irR[\n(8R/a) - f]. Therefore, we see that A depends only on the aspect ratio R/a of the ring, and not on its 
circumference: 



A « 6.0x10 



-5 



ln(^ 

a 



(42) 



For a thick ring, where R/a « 1, C is approximately equal to the size L of the ring. For a thin ring with R/a « 10 
or i?/a w 100, £ is substantially larger. 

Consider, for example, a gold ring with L ss 12/im and R/a w 30, characteristic of the rings studied by Chan- 



dresakhar ei a/, [j4|, where persistent currents of order I^ a . 
negligible reduction that is consistent with the experiments. 



where measured. Then C/L « 7.5, and A 
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FIG. 1. Energy difference between the ground state if) and the rotating state ip', as a function of the imparted angular 
momentum 8L. 
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